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The ground-state wave function and the energy gap are calculated for various layer separations, 
d, and for up to 24 electrons by the density matrix renormalization group (DMRG) method. 
Two-particle distribution function and excitonic correlation function are calculated from the 
ground-state wave function, and the evolution of the ground state with increasing d is analyzed. 
The results indicate that the transition is continuous. A smooth crossover of the ground state 
is found at d/l ~ 1.6 from the excitonic character at small d/l to the independent Fermi-liquid 
character at large d/l, where I is the magnetic length. 
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Two-dimensional electron systems under high mag- 
netic field have many interesting features. This is typ- 
ically shown in the ground-state phase diagram 1 ' of 
monolayer quantum Hall systems: incompressible liq- 
uids, 2 ' compressible liquids, 3, 4 ' CDW states 5 ' called 
stripes, bubbles and Wigner crystals are realized only 
by applying uniform magnetic field. 

In bilayer quantum Hall systems, much more interest- 
ing properties are expected owing to the additional de- 
grees of freedom of layers. 6 ' The excitonic phase, namely 
Haplerin's ^1,1,1 state, is one of the ground states real- 
ized in bilayer quantum Hall systems at the total filling 
v = 1 at a small layer separation, d, where electrons and 
holes in different layers are bound to each other due to 
strong inter layer Coulomb interaction. 7,8 ' This excitonic 
state has recently attracted attention because a dramatic 
increase in zero bias tunneling conductance between the 
two layers, 9 ' and the vanishing of the Hall counterflow 
resistance are observed. 10,11 ' 

As the layer separation is increased, the excitonic 
phase vanishes, and at a sufficiently large separation, 
the composite-fermion Fermi-liquid state is realized in 
each layer. Several scenarios have been proposed for the 
transition of the ground state as the layer separation in- 
creases, 12 ~ 18 ' but how the excitonic state develops into 
the independent Fermi-liquid state has not been fully un- 
derstood. In this paper we investigate the ground state 
of v = 1 bilayer quantum Hall systems using the density 
matrix renormalization group (DMRG) method. 19 ' We 
calculate the energy gap, two-particle distribution func- 
tion g(r) and excitonic correlation function for various 
values of the layer separation d, and show the evolution 
of the ground state with increasing d. The results indi- 
cate that the transition of the ground state from d = 
to 00 is continuous. There is a smooth crossover of the 
ground state at around d/l ~ 1.6 from the excitonic char- 
acter to the independent Fermi-liquid character. 

The Hamiltonian of the bilayer quantum Hall systems 



is written as 

H = e-^/V^^-i?!,) 
i<j q 

i<j q 

+ EE^«) e-^e-^/y^M-i^), (i) 

ij q 

where Ri,i are the two-dimensional guiding center coor- 
dinates of the ith electron in layer-1 and i?2,j are those 
in layer-2. The guiding center coordinates satisfy the 
commutation relation, [Rj,Rf\ — il 2 Sjk. V{q) = 2ne 2 /q 
is the Fourier transform of the Coulomb interaction. We 
consider uniform positive background charge to cancel 
the component at q = 0. The wave function is projected 
onto the lowest Landau level. We assume zero interlayer 
tunneling and fully spin-polarized ground state. 

We calculate the ground-state wave function by the 
DMRG method, 19 ' which is a real space renormalization 
group method combined with the exact diagonalization 
method. The DMRG method provides the low-energy 
eigenvalues and corresponding eigenvectors of the Hamil- 
tonian within a restricted number of basis states. 20 ~ 22 ' 
The accuracy of the results is systematically controlled 
by the truncation error, which is smaller than 10~ 4 in 
the present calculation. We investigate systems of vari- 
ous size with up to 24 electrons in the unit cell keeping 
500 basis states in each block. 

In Fig. 1 we first show the lowest excitation gap at the 
total filling factor v = 1 . The total filling factor is defined 
by v — N e /N, where N e — Ni + N2 is the total number 
of electrons and N is the number of one-particle states 
in each layer. N\ and N2 are the numbers of electrons 
in layer-1 and layer-2, respectively, and we impose the 
symmetric condition of Ni = N2. The gap shown in 
Fig. 1 is calculated for N e = 24 and L x /L y — 1.6 in the 
unit cell L x x L y with the periodic boundary conditions 
in both x and y directions. The aspect ratio L x /L y is 
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Fig. 1. Lowest excitation gap A of bilayer quantum Hall systems 
at total filling factor v = 1. N e = 24 and L x /L v = 1.6. The inset 
shows the result for N e = 18 and L x /L y = 1.0. The excitation 
gap A is in units of e 2 /l. 



chosen from the minimum of the ground-state energy 
with respect to L x /L y at around d/l — 1.8, where the 
minimum structure appears in the ground-state energy. 

We can clearly see a large excitation gap for d/l < 1.2, 
where the excitonic ground state is expected both the- 
oretically and experimentally. 7 ~ 18 ) The excitation gap 
rapidly decreases with increasing d/l from 1.2, and it 
becomes very small for d/l > 1.6. This behavior is con- 
sistent with experiment. 23 ^ Although the excitation gap 
for d/l > 1.7 is not presented in the figure for N e = 24 
because of the difficulty in calculating energy of excited 
states in large systems, we do not find any sign of level 
crossing in the ground state up to d/l ~ 4, where two lay- 
ers are almost independent. These results suggest that 
the excitonic state at small d continuously crosses over 
to a compressible state at large d. In the present calcu- 
lation, it is difficult to conclude whether the gap closes 
at a finite d/l ~ 1.6 in the thermodynamic limit or the 
excitonic state survives with an exponentially small fi- 
nite gap even for a large d/l > 1.6. We have calculated 
the excitation gap in systems of different size and aspect 
ratio, and obtained similar results as shown in the inset 
of Fig. 1. 

In the first excited state, however, a level crossing is ex- 
pected at d/l ~ 1.2, where we can see a sudden decrease 
in excitation gap, which is due to a crossing of excited 
states. This means the character of the low-energy ex- 
citations changes at d/l ~ 1.2, which is confirmed by a 
clear change in the correlation functions in the excited 
state as shown later. 

To study how the excitonic ground state at d = 
crosses over to the compressible state at a large d, we 
next calculate the exciton correlation defined by 

ScxM = (*l4,n c 2,n4,0Cl,0|*), (2) 

where \^) is the ground state and c\ n {c\ n ) is the cre- 
ation operator of the electrons in the nth one-particle 
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Fig. 2. Exciton correlation of bilayer quantum Hall systems at 
v = 1. The solid line represents g ex (7V/2). The dashed line 
represents g cx (N/2 — 1). 



state defined by 

Mr) = - / = m exp \ik y y - ^f^} (3) 

in layer-1 (layer-2). X n = nL x /N = —k y l 2 is the x- 
component of the guiding center coordinates and L x is 
the length of the unit cell in the x-direction. 

Because g ox (n) represents the correlations between 
the two excitons created at X = and X — X n , 
liniTj^oo g ex (n) ^ indicates the existence of the macro- 
scopic coherence of excitons. As shown in Fig. 2, g cx (n) 
tends to 1 as d — > 0. Indeed Haplerin's W^i.i state has 
this macroscopic coherence and g cx (n) — 1 independent 
of n. In this figure we have shown g ex (N/2) instead of 
lim, woo g ex (n), because the largest distance between the 
two excitons is L x /2 in a finite unit cell of L x x L y un- 
der the periodic boundary conditions. To check the size 
effect, we also plot g ex (N/2 — 1) with a dashed line. Be- 
cause the difference between g ex (N/2) and g ex (N/2 — 1) 
is small, we expect g ex (N/2) well represents the macro- 
scopic coherence in the limit of N — > oo. 

With increasing d, the excitonic correlation decreases 
monotonically and finally falls down to a negligible value 
at d/l ~ 1.6. Since g ax (N/2) continuously approaches 
zero, the transition is not of the first order. This sec- 
ond order or crossover transition at around d/l ~ 1.6 is 
consistent with the result of the excitation gap shown 
in Fig. 1. The slight difference in the asymptotic form 
around d/l ~ 1.6 is expected to be a finite size effect. 

We next calculate pair correlation functions of the elec- 
trons to see the detailed evolution of the ground-state 
wave function. The interlayer pair correlation function 
of the electrons is defined by 

Sl2(r) = ^-{M E + - fl2,m)|*>. (4) 

^ n m 

We present Agi2(r) in Fig. 3, which is denned by 

A 5l2 (r) = / ( 5l2 (r') - l)*(|r'| - r) dr', (5) 
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Fig. 3. Interlayer pair correlation function of electrons in ground 
state of bilayer quantum Hall systems at v = 1. N e = 24 and 
L x /L y = 1.6. Inset shows 2Api2(r)/r, which is essentially the 
same as gi2{r) — 1. 



Fig. 4. Intralayer pair correlation function of electrons in ground 
state of bilayer quantum Hall systems at V = 1. N e = 24 and 
L x /L y = 1.6. The inset shows 2Agn(r)/r, which is essentially 
the same as gn(r) — 1. 



where r' is the two-dimensional position vector in each 
layer. Ag 12 (r) represents the difference from the uniform 
correlation of independent electrons multiplied by 2irr. 
We also show 312 (r) — 1 in the inset. 

At d — we find clear dips in Agi 2 (r) and gi 2 (r) at 
r/l ~ 1 and 0, respectively, which are a characteristic 
feature of the excitonic state made by the binding of 
electrons and holes between the two layers. The binding 
of one hole means the exclusion of one electron caused by 
the strong interlayer Coulomb repulsion. The increase in 
layer separation weakens the Coulomb repulsion between 
the two layers and reduces |Ag 12 (r)| at r/l ~ 1. 

The decrease in the interlayer correlation |Agi2(r) 
opens space to enlarge the correlation hole in the same 
layer and reduce the Coulomb energy between the elec- 
trons within the layer. This is shown in Fig. 4, which 
shows the pair correlation functions of the electrons in 
the same layer defined by 

9n(r) = jvJ^TT) (*l £ S ( r + R ^ ~ «i,m)l*X6) 

nm 

A511M = / Cffu(r') - l)5(W\ - r) dr'. (7) 

The obtained results indeed show that the correlation 
hole in the same layer at r/l — 1 is enhanced with the 
increase in d contrary to the decrease in the size of the 
interlayer correlation hole in Fig. 3. The correlation hole 
in the same layer monotonically increases in size up to 
d/l = 1.8, and then it becomes almost constant. The 
correlation function gn(r) for d/l > 1.8 is almost the 
same as that of v = 1/2 monolayer quantum Hall systems 
realized in the limit of d = 00. This is consistent with the 
almost vanishing excitation gap and exciton correlation 
at d/l > 1.8 shown in Figs. 1 and 2. 

In the v = 1/2 monolayer system, the ground state 
is described as a dipolar composite fermion liquid state, 
where an electron forms a bound state with a displaced 



hole. 24 -' The peak of A<7n(r) at r/l ~ 4 shows the posi- 
tion of the nearest neighbor electrons in accordance with 
the approximate mean distance between the electrons 
3.54Z estimated from {L x L y /Nx) x / 2 = {2ttN/N 1 )^ 2 1. On 
the other hand, the correlation hole at r/l ~ 1.5 shows 
the existence of the bound hole forming the dipolc. 

As the distance between the two layers decreases, the 
interlayer correlation between the dipoles begins to de- 
velop. This is seen in Fig. 3 as the development of the 
correlation hole at r/l ~ 1 and a peak at r/l ~ 2.5 
in A<7i2(r). That Apn(r) remains almost the same at 
d/l > 1.6 indicates that the identity of the dipolar com- 
posite fermion remains in this region. The almost same 
amplitude of Agi 2 (r) at r/l ~ 1 and 2.5 for d/l > 1.6 is 
also consistent with the antiparallel interlayer correlation 
of dipolar composite fcrmions. 

With further decrease in d/l from 1.6, interlayer corre- 
lation between electrons and holes becomes stronger and 
begins to compete with intralayer correlation. The dipo- 
lar composite fermion begins to lose its identity, and the 
ground state gradually changes into a collection of in- 
terlayer electron-hole pairs, i.e., interlayer excitons. The 
dissociation of the dipole is seen in the rapid decrease in 
the size of the intralayer correlation hole at d/l < 1.2 in 
Fig. 4, and in the broadening of the peak at r/l ~ 2.5 in 
Fig. 3. On the other hand, the development of the cor- 
relation hole at r/l ^ 1 in Fig. 3 and the development of 
the exciton correlation in Fig. 2 shows the formation of 
interlayer excitons. The dissociation of the dipole is also 
considered as continuous weakening of the antiparallel 
correlation between interlayer excitons as d decreases. 
This correlation disappears at d = when the ^1,1,1 
state is realized, and Apn(r) and Agi 2 (r) coincide. In 
this course, half of the correlation hole in the same layer, 
which is a characteristic feature of composite fermions, 
is transferred into the other layer. 

The breakdown of the composite fermions at d/l ~ 1.2 
affects the character of the lowest excitations, which is 
clearly shown in the level crossing of the excited state at 
d/l ~ 1.2. The change in the character of the excitation 



Naokazu Shibata and Daijiro YOSHIOKA 



0.01 



0.005 



-0.005 



-0.01 



---- °£n 
" 5g 12 






d/l=0 2 


0.4 

\_0.6 


1.4 

- #Gw' 
/ ** \\ - 


13 ^ v 






V\.0.8 

XVV 1 - 
1.2 






// 




,--;,^-i.o - 

-''7^-0.8 










7^0.6 








dlf=0.2 


-0.4 



4 

r/Z 



Fig. 5. Change in correlation function through excitation from 
ground state to first excited state, v = 1 and N e = 18 with 
L x /L y = 1.0. 



is confirmed by the correlation functions in the excited 
state. Figure 5 shows the difference in the pair correla- 
tion functions gij{r) between the ground state and first 
excited state defined by 



L X Ly 



(9?Ar')-9?Ar'))5(\r'\-r) dr',(8) 



where gfj{r) and gfj(r) are the pair correlation functions 
in the ground state and the first excited state, respec- 
tively. 5gij (r) in Fig. 5 shows that there is a discontinu- 
ous transition between d/l = 1.2 and 1.3, which supports 
the level crossing in the first excited state. 

Below d/l ~ 1.2, Sgij(r) have a large amplitude at 
r/l ~ 2 and 6. In this region of d/l, interlayer-excitons 
are formed and they are weakly interacting. Figure 5 
shows that the antiparallel correlation between excitons 
realized in the ground state is weakened in the excited 
state. A small singularity at r/l ~ 5.3 is due to the finite 
size effects of the square unit cell. Above d/l ~ 1.2, only 
^ffi2(?") have large amplitude at r/l ~ 1 and 2. In this 
region, dipolar composite fermions are formed and there 
is interlayer correlation among them. Figure 5 shows 
that the low-energy excitations are made by modifying 
this interlayer correlation between dipoles. 

In summary, we have investigated the ground state and 
low-energy excitations of v = 1 bilayer quantum Hall 
systems by the DMRG method. The obtained exciton 
correlations show that the ground state for small layer 
separation, d/l < 1.6, is characterized by the macro- 
scopic coherence of excitons, g cx ^ 0, which is consis- 
tent with the vanishing of the Hall countcrflow resis- 
tance 10 ' 11 ^ and a marked enhancement in zero-bias tun- 
neling conductance 9 ^ at low temperatures. The excita- 
tion gap A rapidly decreases with increasing d/l from 1.2 
to 1.6, which is also consistent with the clear decrease in 
zero bias tunneling conductance 9 ' and transport experi- 
ments 10, n > 23 ) at finite temperatures around this region. 

We have also studied the evolution of the ground-state 
wave function. We have found a clear dip in A<?i2(r) 
at r/l ~ 1 and a peak at r/l ~ 2.5 for d/l > 1.2, 
which is consistent with the interacting dipolar compos- 



ite fermion picture. 15 - ) At d/l ~ 1.2, the peak in Agi2(r) 
at r/l ~ 2.5 becomes largest, and the correlation func- 
tions becomes almost the same as those of the attractive 
Vj 12 model. 17 - ) The effective interlayer attractive interac- 
tion between the electrons is induced by the antiparallel 
correlation between dipolar composite fermions, which 
can also be interpreted as the antiparallel correlation be- 
tween interlayer excitons. A further decrease in d con- 
tinuously transfers the correlation hole to the other layer 
and the interacting interlayerexcitons gradually become 
a better picture of the ground state. We do not find any 
sign of the first-order transition through the above evo- 
lution of the wave function. Thus, the transition from 
the compressible state at d = oo to the excitonic state 
at d = is continuous. 
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